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Abstract 

Based on a Glauber model, a statistical analysis of all mid-rapidity J/ip 
hadroproduction and leptoproduction data on nuclear targets is carried out. This 
allows us to determine the J/'i/'-nucleon inelastic cross section, whose knowledge 
is crucial to interpret the J/ijj suppression observed in heavy-ion collisions, at 
SPS and at RHIC. The values of cTj^^^ are extracted from each experiment. A 
clear tension between the different data sets is reported. The global fit of all data 
gives Cj/^,i>j = 3.4ib0.2 mb, which is significantly smaller than previous estimates. 
A similar value, o""^^ = 3.5 it 0.2 mb, is obtained when the nDS nuclear parton 
densities are included in the analysis, although we emphasize that the present 
uncertainties on gluon (anti) shadowing do not allow for a precise determination 
of fj^^jj. Finally, no significant energy dependence of the J/ip-N interaction is 
observed, unless strong nuclear modifications of the parton densities are assumed. 



1 Introduction 



Twenty years ago, Matsui and Satz suggested the suppression of heavy-quark bound 
states in heavy-ion colhsions as a sensitive probe for quark-gluon plasma (QGP) forma- 
tion [1]. Since then, an intense experimental activity has been carried out. Soon after 

§Email address: arleoOcern. ch (corresponding author) 

^On leave from Laboratoire d'Annecy-le-Vieux de Physique Theorique (LAPTH), UMR 5108 du 
CNRS associee a TUniversite de Savoie, B.P. 110, 74941 Annecy-le-Vieux Cedex, France 
II Email address: vntram@lbl.gov 
**Present address: Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley CA, USA 
94720 



1 



the Matsui and Satz proposal, the first measurements of J/V' production in heavy- 
ion reactions were performed in 0-U and S-U colhsions by the NA38 fixed-target 
experiment at the CERN Super Proton Synchrotron (SPS) [2]. These were followed 
a few years later by the more precise NA50 data in Pb-Pb collisions [3], and more 
recently by the NA60 preliminary results in In-In collisions [4], at a similar energy 
(V-^NN — 20 GeV). Finally, the PHENIX collaboration reported on data on J/ if) pro- 
duction in Au-Au [5] and on preliminary measurements in Cu-Cu scattering [6] at 
the Relativistic Heavy Ion Collider (RHIC) energy, ^/s^ = 200 GeV. Remarkably, all 
these experimental results indicate a significant J/ip suppression in heavy ions (prop- 
erly normalized by the number of NN binary collisions) with respect to p-p scattering. 

However, it rapidly became clear that a key requirement to interpret these data is 
the quantitative understanding of all nuclear effects - i.e. those that are not due to QGP 
formation - that could also affect J/ip production in nucleus-nucleus reactions. Among 
them, the nuclear modifications of the parton densities, measured in deep inelastic 
scattering (DIS) and Drell-Yan (DY) reactions, could affect the nuclear dependence of 
J/ ip production. Another effect is the inelastic rescattering of the J/ ip state - or more 
generally a cc pair^ - in cold nuclear matter, which is expected to give the dominant 
contribution to the J/ ip suppression in proton-nucleus and peripheral nucleus-nucleus 
collisions. The strength of this mechanism is monitored by essentially one physical 
parameter, the J/^/^-nucleon inelastic cross section, (Tj^^^- 

On the theoretical side, this cross section has been computed perturbatively in QCD 
in the leading-twist approximation [7] and proves of the order of a few millibarns in the 
J/ip-N energy range of interest [8,9] (similar estimates based on a generalized vector 
dominance model are given in [10]). However, the relative smallness of the charm quark 
mass (with respect to A^^^) may lead to significant power corrections; consequently 
these numbers should not be taken at face value but rather seen as qualitative estimates 
only. It is therefore necessary to rely on existing experimental data to determine this 
quantity. Fortunately, the measurement of J/ ip production on nuclear targets has been 
performed extensively for hadropro duct ion (E537 [11], NA3 [12], NA38 [13], NA50 [14], 
NA60 [15], E672 [16], E772 [17], E866 [18], HERA-B [19], PHENIX [20]), photopro- 
duction (SLAG [21], E691 [22]) and leptoproduction (EMG [23], NMG [24]). In the 
spirit of the Gerschel and Hiifncr analysis [25] and the subsequent work by Kharzeev, 
Lourenco, Nardi and Satz [26], it is the aim of this paper to extract the inelastic J/ip- 
N cross section systematically from a Glauber analysis of these experimental results 
- including the possible nuclear corrections to the parton distribution functions - and 
thus to provide a baseline of the J/ ip absorption mechanism in nuclear matter from all 
available data. 

The outline of the paper is the following. The J/ijj production cross section in 
nuclei is detailed in Section 2. The method adopted in this analysis is then described 

^Throughout this paper, we will abusively mention the "J/^/^-N interaction"; yet at high energy 
what propagates through the nucleus is a cc pair. Similarly, we later denote by "0"^^^,^" the parameter 
governing the suppression of J/ip states in nuclear collisions, whether or not the J/ip is formed inside 
or outside the nuclear medium. The issue on the finite J/tp formation time is addressed in Sect. 7. 
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in Section 3, while Section 4 is devoted to the data selection. The results obtained 
are presented in Section 5. Section 6 is devoted to a systematic investigation of the 
uncertainties coming from the nuclear parton densities, before the results are discussed 
in Sect. 7. Finally, a summary of what has been carried out here is given (Section 8). 



2 Jl^jj production in nuclear targets 
2.1 Production 

The charmonium production cross section in hadron [h) - nucleus (A) collisions is 
determined in this work within the Colour Evaporation Model (CEM). Neglecting the 
possible inelastic interaction of the produced J /ip in nuclear matter, the leading order 
(LO) production cross section, as a function of the longitudinal momentum fraction, 
Xp, reads [27] 



d (jP™'^ 



{hA^ J/i) X) = p 



J/V 



dm 



2m 



■\/x^s + 4m^ 



+ Yl U^^i'^') fti^.,rn') + f^{x„m') f^{x,,m')\ag,im') 



=u,d,s 



where ^ 

^1,2 = 2 +4 "^VS ± ^f) 

are the projectile and target-parton momentum-fractions {^/s being the centre-of-mass 
energy of the hadronic collision), m^ = 1.2 GeV (respectively, = 1.87 GeV) is 
the charm-quark (D-meson) mass, and a^g (respectively, agg) is the LO cc partonic 
production cross section in the quark-antiquark annihilation qq (gluon fusion, gg) 
channel. In this model, the non-perturbative transition from the cc pair to the J/ip 
state is hidden in the proportionality factor, Pj^^ ~ 0.5, which is adjusted to the 
experimental measurements [28]. Note that, since only cross-section ratios of heavy 
to light nuclei will be considered in the following, this factor proves irrelevant in the 
present analysis. In Eq. (1), //' and ff^ stand respectively for the parton distribution 
function (PDF) in the hadron and in the nucleus. The PDF in a nucleus with Z 
protons and the atomic mass number A is written as the sum of the proton (ff^^) and 
the neutron (f^^^) parton densities in a nucleus: 

f^''=z f!"- +{A-z) (2) 

where f^^^ is obtained from f^^^ by isospin conjugation: u^/^ — dP/^^ d"/^ = u^/^^ 
gu/A _ gp/A rpj^-g jgQgpjjj effect only affects the J/t/j production cross section in nuclei 
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marginally, as long as gluon fusion is the dominant channel. However, whenever the 
projectile h contains valence antiquarks (e.g. p, tt^*"), the (f'q^ annihilation process be- 
comes dominant, provided ^/s is not too large with respect to rrij^^ (say, \/s < 20 GeV), 
and significant effects are expected because of the very different valence-quark distri- 
bution in the proton and in the neutron. 

In virtual photon reactions, the J/%Ij inelastic production cross section is computed 
in QCD to leading order according to Refs. [29] : 

d rrP''°^ r 

^(7* A ^ W X) = -1^ /(X, s/m]^^) <(x, Q\ (3) 

with s being the square of the 7*-N centre-of-mass energy, and where Tee — 4.7 keV is 
the J/'^ leptonic width and f{x, s/w?^^) is a sharply peaked function above rri^^^^/s [29]. 
Just as Pj^^ in (1), these two factors cancel out in the heavy-to-light nucleus cross- 
section ratio. We shall therefore only exploit the fact that the J/ip inelastic leptopro- 
duction in nuclei is proportional to the gluon density in the nuclear target, /^(a;, Q^), 
where the hard scale Q is given by the J/il) mass. 

The PDFs in Eqs. (1) and (3) are taken from the LO paramctrizations CTEQ6L 
for the proton [30] and SMRS for the pion [31]. The nuclear PDFs are discussed in the 
next section. 



2.2 Nuclear part on densities 

The nuclear parton densities arc known to differ somehow from those in a proton, 
fi^^i^) 7^ fii^) (^66 [32], and references therein, for a recent review) on the 
whole Bjorken-x range. At small x < 0.05, nuclear PDFs (nPDF) are depleted, 
fi^^{x) < fi{x), in the so-called shadowing region, while they may be slightly en- 
hanced ( "antishadowing" ) in an intermediate x ~ 0.05-0.2 range. Above x > 0.2, the 
quark distributions were determined from the EMC muoproduction data on nuclei [33] 
and proved smaller than the quark densities in a proton (depletion usually referred to 
as the "EMC effect"). At very large a; ~ 1, f^^^{x) :§> ff{x) because of the Fermi 
motion of the nucleons in the nuclei. 

Using DIS and DY measurements on various nuclear targets, several LO fits of 
the nuclear parton densities, analogous to the standard QCD fits in a proton, have 
been performed, first by Eskola et al. (EKS) [34], and more recently by Hirai et al. 
(HKM) [35] and Dc Florian and Sassot (nDS) [36]. Note that in the latter analysis, a 
parametrization at NLO accuracy is also given. Unfortunately, the current precision 
(and variety) of the nuclear data does not yet allow the nPDFs to be as constrained as 
those in the proton, leading in particular to rather large uncertainties (say, 20% in the 
X range of interest here) in the gluon sector. This is illustrated in Fig. 1, where the ratio 
of the gold nucleus over the proton PDF at the J/ip mass scale, R^^{x, — 171^^^^), is 
plotted as a function of x using the different paramctrizations. 



4 



~1 1 1 — I — I — I I I 



"T 1 1 1 1 — I I I 



1.2 



— nDS 
■■■ nDSg 

- EKS 
-■- HKM 




1 1 — r-] — I — r 

/ 

/ 

/ 

/ 



I" 
I 



\ / 



0.6 



0.001 



0.01 



0.1 



Figure 1: The ratio of the gluon distribution in a gold nucleus over that in a proton, 
fn"]^^), plotted as a function of Bjorken-x using the nDS, nDSg, EKS and HKM 
parametrizations. 



Since the nuclear dependence of the gluon and the quark densities are a priori 
not identical, ) 7^ R^[x,rn? ), it is necessary to determine the relative 

proportion of the gg fusion and qq channel in the J hadroproduction process, which 
is given here by the CEM, Eq. (1). However, we believe that this proportion - and 
hence our predictions - should not depend on the specific model assumed for the heavy- 
quarkonium production mechanism (Colour Singlet Model, Colour Evaporation Model, 
Non-Relativistic QCD, . . . ); see e.g. [37]. 

The parton distributions in nuclei not being precisely fixed yet by the current 
experiments, the analysis is first carried out without any nuclear corrections to the 
PDFs, assuming fi^^{x) = fj{x). The analysis is also performed using the LO nDS 
parametrization, whose agreement with DIS and DY nuclear data is the best among 
the different nPDF sets [36]. Furthermore, we shall critically discuss in Section 6 how 
the uncertainties on gluon ( ant i) shadowing affect our results. 



2.3 Nuclear absorption 



The factorization between the charmonium production process and the subsequent 
possible J/ip inelastic interaction with nuclear matter is assumed. Quite generally, we 
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may thus write the J/%1) production cross section as 

^ % 7* A ^ J/V^ X) = 5,bs(A , a,/, J X {h, 7* A ^ J/^^ X) , (4) 

where 5'abs(A ,cr^^^j^) denotes the probabihty for no interaction (or "survival proba- 
bihty") of the J/V' meson with the target nucleus. It depends of course on both the 
atomic mass number A of the nucleus and the J/ip-N inelastic cross section, cTj^^^, and 
is given in a Glauber model by [38] 

S.^(A..„„)= (^3^^/ db(l-e-'-/-'"-'^<-"^«"), (5) 
with the thickness function T^(b) 

/+0O 
dzp{h,z). (6) 
-oo 

The thickness function Eq. (6) is determined for all nuclei using a 2-parameter Fermi 
model for the nuclear density profile p{h,z), except for the Be, C, and Ca targets 
where a harmonic oscillator, a sum of gaussians, and a 3-paramctcr Fermi model is 
used respectively (see Ref. [39] for details) and for the Ag and Ft nuclei whose density 
profiles, missing in the analysis [39], are taken from the parametrization given by the 
FRITIOF Monte Carlo code [40]. Finally, we neglect the nuclear absorption in the 
proton and deuterium targets, that is 5'abs(^ < 2 , cTj. ) — 1. 



3 Method 

After discussing the theoretical calculations of J/ip hadroproduction and leptoproduc- 
tion on nuclei, we present in this section the method followed in this study. The J/ip-'H 
inelastic cross section is extracted, for each experimental sample £ with data points, 
from the minimization of the function [41] 

2 

- VM-^V, (7) 
where the theoretical nuclear production ratio, 

(/i,7*B^ J/V'X) (8) 

X. 

% 

computed from Eq. (4) at a given x^. or Bjorken a;., depends explicitly on the free, 
but positive, parameter, o ^^^^ (which we denote generically by c^^'^^ when nuclear 
parton distributions are used in the theoretical calculations). In Eq. (7), (7j represents 
the statistical and the uncorrelated systematic errors, added in quadrature, of the 



B da 
A dx 



(/i,7* K^Jli= X) 



do; 
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experimental nuclear production ratio, Rf^"^, in a heavy nucleus (A.) with respect 
to the lightest target available in the data sample^ (B). In order to avoid too large 
systematic errors in the experimental ratio, both reactions on targets A. and B are 
required to be taken from the same experiment and at the same centre-of-mass energy. 

The correlated systematic error from source k on the data point i, Pik, appears in 
(7) through the i^-dimensional vector 



and the K x K matrix 



Mki = 5ki + Y. (10) 

where K is the number of distinct correlated errors. In the present work, we shall 
mostly be concerned by normalization errors (e.g. due to the uncertainty on the lightest 
nucleus production cross section), and the number of distinct errors in eac/i experiment 
is K = 1. In practice, M is a number when (y j^^^ is extracted from one experiment, 
and a diagonal M x M matrix when performing a global fit on the M data samples. 

The definition, Eq. (7), is shown to be actually equivalent (see [41]) to the more 
traditional approach, which involves the inversion of the n x n covariant matrix: 

K 

W = + E (^^^^ P^k- (11) 

fc=l 

Since here K <^ n, we will use the former definition (7) in the present paper, which 
proves of course more efficient than Eq. (11). The phenomenological consequences of 
taking the correlated errors properly into account in this analysis - instead of adding 
all errors quadratically - will be discussed in Sect. 5. 

As already stressed, the estimated J/ip-N cross section, d-j^^^, is determined from 
the minimization of the function: 

x'('^jmn) = min x'(a^,^^) = . (12) 



Including the correlated errors in the definition, the la error 5a j^^^ on the fitted 
parameter^ ^j/v,n l^ads to a deviation of by one unit from its minimum: 

Ax' = x'(%v,N ± ^^j/*n) - xL = 1- (13) 

Note, however, that this criterion is no longer valid when correlations are neglected; in 
which case, it should be emphasized that Ax^ can actually be much larger than 1 [41]. 



^Notc that for the NA3 experiment, the inverted ratio I^l^'^ (p /Vi) is measured. 

■^We drop the subscript "min" in the following, and whenever a value is quoted, it is understood 
to be xLn- "^^^ "hat" on the estimated cross section aj^^^, and its error S&j^^j^, will also be removed 
in the following, for clarity. 
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For some data samples, the agreement between the theoretical predictions and the 
experimental measurements is poor, leading to rather large ^ 1 values. Following 
the prescription from the Particle Data Group [42], the la error is rescaled. 




^ S X 5a 



(14) 



where the factor S is defined as: 



11/2 



s = 



if xVndf>l, 



(15) 



and S = 1 otherwise. 

The J/ijj-'N cross section is systematically determined in Sect. 5 from the individual 
data samples. It will also be extracted from a global fit of all available data, from the 
minimization of the weighted function: 



This global fit analysis will thus give lesser importance to the data sets for which no 
good agreement between data and theory is observed. 

4 Data selection 

The data sets used in this analysis are summarized in Table 1, in which the centre-of- 
mass energy ^/s of the hadron-nuclcon or photon-nucleon system, the x^, domain^, the 
projectile/target species and the number n of points in each data sample are specified. 

4.1 Hadroproduction data 

Concerning hadroproduction measurements, we analyse the data from the E537, NA3. 
NA38, NA50, E672, E866, HERA-B, and PHENIX collaborations. The projectiles 
used were mainly protons (NA3, NA38, NA50, E866, HERA-B, PHENIX), but also 
pions (E537, NA3, E672), antiprotons (E537), and deuterium nuclei (PHENIX). The 
centre-of-mass energy range is ^/s — 15-200 GeV which, in the x^, domain covered by 
each experiment and which fits the kinematic requirements detailed below, translates 
into an x,^ = 0.015-0.3 window for 97% of the data points^. Note that the energy in 
the J/V^-nucleon system, given by 



^Bjorken x is understood for the leptoproduction experiment. 

^i.e. excluding the two forward measurements by PHENIX at a;^ = 2 10~^ and = 3.6 10~^. 




(16) 




(17) 
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Exp. 


Ref. 


(GeV) 




Proj. 


Target 


No. data 


E537 


[11] 


15.3 


0-0.25 


TT", p 


Be, Cu, W 


9 


NA3 


[12] 


16.7-22.9 


0-0.3 


7r±, p 


p, Pt 


21 


NA38 


[13] 


29.1 


0.1 


p 


C, Al, Cu, W 


3 


NA50 


[14] 


27.4-29.1 


-0.075-0.075 


p 


Be, Al, Cu, Ag, W, Pb 


13 


E672 


[16] 


31.5 


0.3 




C, Al, Cu, Pb 


3 


E866 


[18] 


38.1 


-0.1-0.3 


p 


Be, Fe, W 


10 


HERA-B 


[19] 


41.6 


-0.3-0.1 


P 


C, W 


6 


PHENIX 


[20] 


200 


-0.1-0.1 


p, d 


p, Au 


5 


NMC 


[24] 


8-21 


0.05-0.15 


7* 


C, Sn 


6 



Table 1: Data sets selected for the present analysis. 



has a rather narrow coverage, ^/sJJ^ — 6.5—25 GeV, despite the wide range in ^/s 
probed by the various experiments. 

The E772 data [17] were not taken into account in this analysis: the incomplete 
coverage of the E772 spectrometer apparatus induces, at small Xp, a strong correlation 
between the longitudinal (say, x^) and the transverse (p^^) momentum of the produced 
J/ip^s. Consequently, J/ip production occurs in this domain mainly at very small Pj_, 
for which the nuclear dependence is the strongest because of the Cronin effect ( [43] and 
M. Leitch, private communication). This explains why the uncorrected J/ip suppression 
reported by E772 in [17] is systematically more pronounced than the more precise and 
more recent E866 data, in the same kinematic domain, which are corrected for this 
bias [18]. 

We would like to mention that the NA60 experiment has recently reported on a 
preliminary measurement of J/ip production in p-A collisions at y/s^ = 17.3 GeV [15], 
i.e. the very energy of the Pb-Pb collisions in which an anomalous J/ip suppression has 
been measured by NA50 (see Introduction). Unfortunately, no J/ip production ratio 
is available yet, which does not allow the present NA60 data point to be used in 
the present analysis. Obviously, the coming measurements in Be, Al, Cu, In, W, Pb, 
and U targets will shed important light on the nuclear dependence of J/i/j suppression, 
and should be used in a future similar analysis. 

Let us finally comment on the kinematic requirements imposed on the available 
hadroproduction data sets. The strong J/ ip suppression reported at large Xp in various 
hadroproduction experiments [12, 17, 18] is so far poorly understood. It may come from 
an interplay of several mechanisms such as inelastic interaction with nuclear matter, 
shadowing effects, and the onset of a new QCD production mechanism coming from 
the intrinsic charm-quark content of the projectile hadrons (see [44], for instance, for 
a quantitative analysis of these various nuclear effects). In any case, it is unlikely 
that inelastic interaction proves the dominant process responsible for such a strong 
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J/ip suppression. Furthermore, since the main motivation of our study is to have a 
handle on the expected amount of J/ip nuclear absorption in heavy-ion collisions, in 
which J/ip^s are mainly produced at mid-rapidity, the large- region is discarded. 
Consequently, an arbitrary cut < 0.3 is therefore imposed on all hadroproduction 
data sets. 

4.2 Leptoproduction data 

As for leptoproduction experiments, the NMC data [24] are analysed in this study. 
Using a 200 GeV and 280 GeV incident muon beam, the virtual photon energy u 
ranges from 40 to 240 GcV in the laboratory frame, corresponding roughly to 7*-N 
centre-of-mass energies ^/s = 8-21 GeV. The typical Bjorken-x range probed in the 
gluon distributions of the nuclear targets is a; = 0.05-0.15 in the NMC kinematics. In 
order to remove the contribution of the J/ ip coherent production, in which the nucleus 
as a whole interacts with the virtual photon, a minimal transverse momentum cut, 
> 0.4 GeV^, was applied. In addition, a cut on the energy fraction carried by 
the J/ip in the laboratory frame, z = Ej^^/u < 0.85, was applied by NMC to safely 
exclude quasi-clastic events [24]. On the contrary, the measurements from the EMC 
collaboration [23] are not taken into account in this work, since their resolution in z 
and did not allow them to subtract the coherent contribution and thus extract the 
genuine J/ip inelastic production, Eq. (3). 

4.3 Photoproduction data 

Regarding the photoproduction experiments, only quasi-elastic events were measured 
at SLAG [21], which prevents us from using these data since the focus is put on J/t/j 
inelastic production. Similarly, the lack of a z energy-fraction cut in the E691 ex- 
periment at Fcrmilab docs not allow for the extraction of the nuclear dependence of 
the inelastic events. Therefore, these measurements are not selected in the present 
analysis. Understanding quantitatively the nuclear dependence of J/ip quasi-elastic 
photoproduction (and leptoproduction) would certainly be interesting in its own, and 
will discussed elsewhere [45] . 



5 Results 

5.1 Determination of (7^/^^ from each experiment 

The results obtained following the analysis described in Sect. 3 are summarized 
in Table 2. The values of the J/ip-N cross section extracted from the individual fits 
of each data sample, as well as their corresponding (rescaled) error, per number of 
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degrees of freedom (x^/ndf) and S factors, are listed in Table 2, using either the proton 
parton density (left) or the nDS nuclear parton density (right). The results using the 
proton PDF are also plotted in Fig. 2, in which the dashed band indicates the J/t/^-N 
cross section determined from the global fit of all data samples. 



Exp. 


(mb) 


xVndf 


S 


J/ipN ^ ' 


„/ndf 

-^nDS' 


S 


E537 


6.6 ± 1.1 


2.0 


1.4 


7.0 ± 1.2 


2.1 


1.5 


NA3 


2.7 ± 0.2 


1.5 


1.2 


2.9 ± 0.2 


1.3 


1.2 


NA38 


5.5 ± 0.7 


3.4 


1.8 


5.4 ± 0.7 


3.5 


1.9 


NA50 


4.7 ± 0.5 


0.4 




4.6 ± 0.5 


0.4 




E672 


11.2 ± 6.1 


0.6 




10.3 ± 5.8 


0.7 




E866 


2.7 ± 0.9 


2.3 


1.5 


2.4 ± 0.7 


1.6 


1.3 


HERA-B 


2.1 ± 1.3 


0.3 




2.0 ± 1.3 


0.2 




PHENIX 


3.5 ± 3.0 


1.7 


1.3 


3.1 ± 2.6 


1.4 


1.2 


NMC 


< 0.9 


0.7 




< 0.8 


0.8 




Global fit 


3.4 ± 0.2 


1.5 


1.2 


3.5 ± 0.2 


1.4 


1.2 



Table 2: The J/t/'-N cross sections, x^/iidf and S factors extracted from each data 
sample with (right) and without (left) nDS nuclear corrections. The bottom row cor- 
responds to the global fit of all the data. 

The results obtained without any nuclear modifications in the PDF are fist dis- 
cussed. From Fig. 2, it is clear, first of all, that not all experiments equally constrain 
the J/V^-N cross section. We may remark in particular that the NA3 and NA50 ex- 
periments allow for a precise determination of cr^^^j^, while the E672, PHENIX and 
preliminary HERA-B measurements lead to a rather large uncertainty in the param- 
eter estimation. It is also worth while to comment on the agreement between the 
theoretical calculations and each data set, which is good (with a low x^/ndf) in the 
E672, NA50, HERA-B and NMC data analysis, but rather poor in the case of the E866 
and, especially, NA38 experiments (x^/ndf = 2.3 and x^/ndf = 3.4, respectively); see 
Table 2. In order to give the reader a feeling on the constraints brought by each 
experiment and on the agreement between data and theory, the relative uncertainty, 
plotted versus the x^/ndf value obtained for the various data sam- 
ples (Fig. 3). We notice, in particular, that the 7% and 11% relative uncertainty in 
the NA3 and the NA50 data is very good, even though x^/ndf remains not too large 
(this is also true, to a lesser extent, with the E537 data). This is in sharp contrast 
with the E672 and HERA-B experiments, whose agreement with theory is good but 
which do not constrain cTj/^^-, and, conversely, with the NA38 experiment for which 
^^j/*n/%*n - 7% and xV^df is very large. 

We now comment on the values extracted from the various data sets. Surprisingly, 
one of the most precise determinations of the J/ V'-N cross section is due to the NA3 
measurements, the "oldest" data set analysed here, thanks to the pretty small system- 
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Figure 2: Values of the J/ip-N inelastic cross section extracted from each data set. 
The solid (dashed) horizontal bars indicate the (rescaled) la error, ScTj^^-^ {6d-j^^^), 
and the dashed band represents the value obtained from the global fit of all data. 



atic and statistical error bars on these data points. As can be seen in Fig. 2, where 
all extracted cross sections are plotted, the NA3 fitted value cTj^^^ = 2.7 ± 0.2 mb is 
perfectly consistent with the estimates from E866 [cTj^^-^ = 2.7 it 0.9 mb), from HERA- 
B {(Jj^^j^ = 2.1 ± 1.3 mb), and with the somehow less precise cr j^^^^ = 3.5 ± 3.0 mb 
from the PHENIX experiment. As already pointed out in [14], a significantly larger 
value, Cj/^N ~ '^■'^ ^ ^-^ ^^^5 determined from the NA50 measurements (this value is 
consistent with (Jj^^^^ = 4.5 ± 0.5 mb extracted by NA50 in [14]). Similarly, a 4-6aj^^^ 
discrepancy is observed between the NA3 central result and the extracted values from 
the E537 and NA38 data (cr^^/^^ ~ 6.6±1.1 mb and cTj/^j^ = 5.5±0.7 mb); yet the rather 
poor data-theory agreement in these samples, possibly due to underestimated system- 
atic errors, weakens somehow this observation (Table 2). The analysis also shows that 
the muoproduction data reported by NMC is basically consistent with no J/ijj nuclear 
absorption, with a not so large la error, Saj^^^ = 0.9 mb and a low x^/ndf = 0.7. 
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Figure 3: The relative uncertainty of the J/ip-'H inelastic cross section, ScTj^^^/aj^^j^, 
plotted versus the values of x^/ndf obtained from each data sample. 



5.2 On correlated systematic errors 



The proper statistical treatment of correlated errors in the present analysis has been 
detailed in Sect. 3. Neglecting such correlations by adding correlated and uncorrelated 
errors in quadrature - yet this is often done in fitting procedures - indeed leads to 
strongly biased values for the J/ip-'N cross section. 

Let us illustrate this with two examples, analysing for instance the E866 and the 
PHENIX data, which both have normalization (i.e. correlated) errors on the ratio 
^exp _ p _ X2%, respectively). The estimate from the E866 experiment 

goes from cTj^^^ = 2.7 ± 0.9 mb (Table 2) to cr^^^j^ = 3.1 ± 0.3 mb when correlations 
are neglected. Similarly, cr j^^^ determined from PHENIX is shifted from (Tj^^^^ = 
3.5 ± 3.0 mb to (Jj^^^^ = 2.0 ± 1.3 mb. The values are thus affected by 15% and 
75% - downwards for E866 and upwards for PHENIX - when correlated errors are 
considered. 

This can be better understood within the "pull approach" , which is another equiv- 
alent way of dealing with correlated systematic errors. When extracting Cj^^^ within 
this alternative approach (with of course identical results), we found that the pulls 
C, (that is the amounts, normalized by (3, by which the theory predictions need to be 
shifted to best accommodate the data) extracted for both experiments indeed are quite 
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large and have opposite signs, ^^ggg = —0.3 and ■Cphenix — Perhaps even more im- 
portantly, the error on the parameter is smaller by a factor of almost 3 when correlated 
errors are not properly taken into account. 

5.3 Global fit analysis 

As discussed in Sect. 5.1, a significant tension between the different data samples is 
thus observed. On the one hand, the values extracted from the NA3, E866, HERA- 
B and PHENIX experiments are somehow consistent with each other (despite some 
spread), while we remark that the E537, NA38, and NA50 results clearly stand some- 
how above the bulk of these estimates. The origin of this discrepancy is not clear to 
us. The vanishing cross section determined from J/ ip muoproduction also stands apart 
from most of the estimates, even though Sa ^^^.^ = 0.9 mb. A significantly different 
nuclear dependence of J/ip hadroproduction and leptoproduction may actually ques- 
tion the factorization assumed in Sect. 2.3. However, more leptoproduction data would 
obviously be needed to clarify this issue. 

Despite the above-mentioned tension reported between the individual data sets, 
the J/V'-N cross section is obtained from a global fit of all data, properly weighted 
to account for the possible discrepancy between each data sample and the theory (see 
Sect. 3). The global fit leads to a combined cross section of cTj^^-^ = 3.4 ± 0.2 mb, and 
naturally a rather large x^/ndf = 1.5. This estimate, as well as its possible energy 
dependence, will be discussed in Sect. 7. 



Excluded expt. 


None 


E537 


NA3 


NA38 


NA50 


xVndf 


3.4 ± 0.2 
1.5 


3.3 ±0.2 
1.4 


4.6 ±0.3 
1.1 


3.2 ±0.2 
1.3 


3.2 ±0.2 
1.5 


Excluded expt. 


E672 


E866 


HERA-B 


PHENIX 


NMC 


(mb) 
xVndf 


3.4 ±0.2 
1.5 


3.4 ±0.2 
1.5 


3.4 ±0.2 
1.6 


3.4 ±0.2 
1.5 


3.4 ±0.2 
1.5 



Table 3: The J/ip-N cross section extracted from the global fit excluding the data from 
one experiment. 

Given that some experiments give stronger constraints on cr^^^^ than others (see 
above and Fig. 3), it is interesting to repeat the global fit excluding each time one 
experiment. The values obtained for cTj^^j^ are quoted in Table 3. It is not surprising to 
sec that the most dramatic changes on (J j^^.y arc observed when the NA38 (despite only 
3 data points), NA50, and essentially the NA3 experiment (from Cj^^^ = 3.4 ± 0.2 mb 
to Uj^^j^ — 4.6 ± 0.3 mb), relaxing significantly the tension (see Fig. 2 and the x^/ndf 
in Table 3). 
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5.4 Using nuclear parton densities 



Let us discuss the results of the analysis when the nDS nuclear parton densities are 
used. We notice in Table 2 (right) that the inclusion of the nPDF in the theoretical 
calculation does not significantly affect the results obtained using the proton PDF. 
Since the shadowing process depletes J/%1) production in nuclei at high energy, the 
extracted inelastic cross sections tend to be slightly smaller, since part of the J/V' 
suppression seen in the data is already accounted for this effect. This is particularly 
true for the highest-energy E866, HERA-B, and PHENIX experiments. We also note 
that the agreement between the E866 and the theory is significantly improved, with 
X^jjg/ndf = 1.6 (and a smaller cross section of c^^^ = 2.4 ± 0.7 mb) with respect 
to x^/ndf = 2.3 {(Tjf^^ = 2.7 ± 0.9 mb) without any nuclear correction in the PDF. 
At RHIC, it is also worth while to remark the better agreement when nuclear parton 
distributions arc assumed in the calculation (x^^^/ndf = 1.4) than when they arc not 
(x^/ndf = 1.7). At lower energy (larger x^), the effect of nDS antishadowing is rather 
tiny (see in particular the increase of cTj^^^ for the E537 and NA3 data samples); 
consequently the values of the extracted cross sections are not dramatically affected 
by these corrections. Note that using nDS PDFs makes x^/ndf closer to 1 for most 
experiments. Fitting all data samples results in a cross section o"^^^^ = 3.5±0.2 mb with 
X^Qg/ndf = 1.4, that is perfectly compatible with cr^^^j^ = 3.4 ± 0.2 mb (x^/ndf = 1.5) 
using the proton densities. 

However, as mentioned in Sect. 2.2, the gluon distributions in nuclei are not tightly 
constrained experimentally (see Fig. 1). It is therefore the aim of the next section to 
discuss further the uncertainties of the nPDFs, and especially the consequences on the 
'^j/V'N values extracted from the data. 



6 Uncertainties from the nuclear parton densities 
6.1 Global fit analysis 

It has ah^cady been stressed that the current DIS and Drell-Yan nuclear measurements 
do not allow stringent constraints to be put on the amount of gluon (anti) shadowing. 
Therefore, in order to quantify the sensitivity of our results with respect to the strength 
of nuclear corrections to the parton densities, the global fit analysis is repeated in this 
section using other nPDFs parametrizations: nDSg [36], EKS [34], and HKM [35]. The 
results are summarized in Table 4. 

The nDSg set is a constrained fit where the gluon shadowing is arbitrarily enhanced, 
setting Rg = 0.75 at = 5 GeV^ in a Au nucleus. As stressed in [36], the agreement 
between the nDSg parametrization and the DIS and DY data is much worse than the 
unconstrained fit nDS (x^^^^^/ndf = 1.4 versus x^^s/^*^^ ~ ^■^)' ^® ^^^^ 
further quantify the systematic uncertainty of the above result, due to the poorly known 
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Proton 


nDS 


nDSg 


EKS 


HKM 


J/rpN ^ ' 


3.4 ± 0.2 


3.5 ± 0.2 


3.9 ± 0.2 


5.2 ± 0.2 


3.5 ± 0.2 


xVndf 


1.5 


1.4 


1.5 


1.6 


1.9 



Table 4: The J/^'-N cross section extracted from the data using the proton and various 
nuclear parton density parametrizations. 



gluon nPDF. We find a somewhat larger cross section, a"^^^^ = 3.9 ±0.2 mb, which is a 
direct consequence of the gluon antishadowing, with a slightly higher x^^^ /ndf = 1.5. 
Even more spectacularly, the stronger antishadowing in the EKS parametrization (see 



Fig. 1) further enhances the J/V'-N cross section up to o ^j^^ — 5.2 ± 0.2 mb. Finally, 
the cross section extracted using the HKM parton densities, cj^^™ = 3.5 ± 0.2 mb, 
is identical to what is obtained with the nDS parametrization; yet the data-theory 
agreement is somehow worse (x^ /ndf = 1.9, as compared with /ndf = 1.4). 



The spread of quoted in Table 4 directly reflects the present lack of knowledge 

of the (gluon) nuclear densities. Clearly, a better determination of o'^fi'^^ could only 
be achieved when these are more tightly constrained by experimental data. However, 
it should be repeated that, among these different parametrizations, the nDS fit offers 
the best agreement with the nuclear DIS and DY data. In that sense, we believe 
"^//^N ~ "^'^ '^■^ most likely estimate. It is also interesting to point 

out that the ^ hierarchy reported here, x^^g < x^sg ~ ^Ls ~ xLm' ^^^^ 
as what is encountered in the DIS and DY context^, although the complexity of the 
Jji) production process (as compared with the "cleaner" DIS and DY probes) makes 
it difficult to draw quantitative conclusions on the reliability of the various nuclear 
distributions. Nevertheless, we shall see in the next section that some experiments 
actually bring stringent constraints on the nPDFs. 



6.2 Systematics of nuclear parton densities 

For completeness, we give in this section the J/V'-N inelastic cross sections (and 
X^/ndf) extracted from the individual data samples using the various fits of the nu- 
clear parton densities. In Table 5 are listed the estimates using the default nDS fit, as 
well as the HKM, nDSg, and EKS parametrizations. These last two exhibit a stronger 
(anti)shadowing than nDS and HKM. 

As already emphasized, it is in the E866, HERA-B and PHENIX kinematical do- 
main that shadowing corrections prove the strongest. Consequently, the main differ- 
ences between the various nuclear distributions occur specifically in this region. The 
results are summarized in Table 5. Overall, the cross sections obtained with HKM 
are rather similar to those of nDS, while EKS estimates are way different. The nDSg 

^De Florian and Sassot found x^^^/ndf = 0.8, x^^^^/ndf = 1.4, Xeks/'^^^ ~ ^'^^ ^ significantly 
larger x^^^^/ndf 1 (Ref. [36], and private communication). 
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nDS 






nDSg 




Exp. 


nDS 




S 


(7"°Sg (mb) 




S 


E537 


7.0 ± 1.2 


2.1 


1.5 


7.9 ± 1.2 


2.2 


1.5 


NA3 


2.9 ± 0.2 


1.3 


1.2 


3.4 ± 0.2 


1.4 


1.2 


NA38 


5.4 ± 0.7 


3.5 


1.9 


6.2 ± 0.8 


3.9 


2.0 


NA50 


4.6 ± 0.5 


0.4 


— 


5.4 ± 0.5 


0.6 


— 


E672 


10.3 ± 5.8 


0.7 


- 


8.8 ± 5.4 


0.7 


- 


E866 


2.4 ± 0.7 


1.6 


1.3 


3.6 ± 1.8 


7.9 


2.8 


HERA-B 


2.0 ± 1.3 


0.2 




2.6 ± 1.4 


0.3 


- 


PHENIX 


3.1 ± 2.6 


1.4 


1.2 


0.6 ± 1.9 


0.8 




NMC 


< 0.8 


0.8 




< 0.9 


0.7 


- 


Global fit 


3.5 ± 0.2 


1.4 


1.2 


3.9 ± 0.2 


1.5 


1.2 






EKS 






HKM 




Exp. 


EKS (-^b) 




S 


HKM (j^b) 


XHKM/ndf 


S 


E537 


8.2 ± 1.1 


1.9 


1.4 


7.1 ± 1.1 


2.0 


1.4 


NA3 


4.6 ± 0.2 


1.2 


1.1 


3.0 ± 0.2 


1.5 


1.2 


NA38 


7.9 ± 0.8 


3.2 


1.8 


5.5 ± 0.7 


3.4 


1.8 


NA50 


6.8 ± 0.5 


0.3 


- 


4.8 ± 0.5 


0.4 


- 


E672 


11.6 ± 6.3 


0.6 


- 


10.3 ± 5.8 


0.6 


- 


E866 


5.3 ± 1.7 


6.5 


2.5 


2.5 ± 0.8 


2.0 


1.4 


HERA-B 


4.2 ± 1.5 


0.9 




2.3 ± 1.3 


0.2 




PHENIX 


1.3 ± 2.0 


0.6 




1.5 ± 2.3 


1.3 


1.2 


NMC 


< 1.6 


0.5 




< 0.9 


0.7 




G1(^1)m1 fit 


5.2 ± 0.2 


l.G 


1.3 


3.5 ± 0.2 


1.9 


1.4 



Table 5: The J/ip-^ inelastic cross section, x^/ndf and S factors extracted from each 
data sample using the nDS, nDSg, EKS, and HKM parametrizations for the nuclear 
PDFs. 



density stands somewhere in between the nDS/HKM and the EKS fitted values. 

It is interesting to note that the cross section extracted at RHIC is basically con- 
sistent with mb (within la) using the nDSg, EKS^, and HKM distributions, while 
a finite (yet with a large 5a j^^^) = 3.1 ± 2.6 mb, is found when the nDS nPDF 

is assumed. Let us also mention that the x^/ndf extracted is much smaller with 
nDSg and EKS (x^^^^/ndf = 0.8 and x^^^/ndf = 0.6) - hence for the distributions 
with the strongest (anti)shadowing - than with nDS and HKM (x^^^^/ndf = 1.4 and 

xLMMf = i.3). 

However, this statement needs to be balanced. Indeed, we can remark that there is 
a good data-theory agreement (say, with x^/ndf ~ 1-2) for all these nuclear PDFs with 
each data sample, with the notable exception of the E866 measurements, which strongly 



^Note that the 0-3 mb range used by R. Vogt for RHIC phenomenology [37] is perfectly consistent 
with our fitted value, a^f^ = 1.3 ± 2.0 mb. 
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disfavour the nDSg (x^^^^/ndf = 7.9) and the EKS (x^^^^/ndf = 6.5) distributions and 
which, in return, definitely support a tiny ( ant i) shadowing a la nDS, with a significantly 
lower x^j^g/ndf = 1.6. Recall that neglecting any nuclear corrections in the PDF yields 
X^/ndf = 2.3 (Table 2), that is an intermediate value between nDS on the one hand, 
and EKS/nDSg on the other hand. Hence, it appears that the E866 J/ip data give 
interesting constraints on the evolution of the gluon density ratio, B^{x,Q'^ ~ m'^^^^), 
in the x = 10~^-10~^ range. 
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Figure 4: The J/ip-N cross section extracted from each data set and in the global fit, 
using the nDS (cr^^^^ ) versus EKS (a^^^ ) nuclear parton densities. 

In order to stress further the uncertainties on Cj^^^ from the unconstrained gluon 
nPDF, is plotted as a function of cr^^^ in Fig. 4. We remark a systematic offset 

of roughly cr^^^ — a""^^ ~ 2 mb for all but one data sample. As mentioned above, 
only the PHENIX experiment leads to values almost twice as large with the nDS than 
with the EKS nPDF. As a consequence, the PHENIX cr^^ stands significantly below 
the other lower-energy experiments (with the exception of NMC), while the PHENIX 
^nDS (^Qgg not depart much from the NA3, HERA-B and E866 estimates. 
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7 Discussion 



We have presented in the previous section the values extracted for the J/^'-N inelas- 
tic cross section from hadropro duct ion and muopro duct ion data on nuclear targets, 
using the simplest geometrical assumptions of the Glauber model. The combined 
cross section has been found to be <Jj/^^ — 3.4 ± 0.2 mb, and a shghtly larger value 
^nDS ^ 3 5 _|_ Q 2 nib was extracted when the nDS nuclear modifications to the PDFs 

J/i/iN 

were taken into account. In this section, we would like first to compare these results 
with previous phenomenological global analyses in the literature, and then with theo- 
retical expectations. The possible energy dependence of the J/V'-nucleon interaction 
is finally addressed in Sect. 7.3. 



7.1 Comparing with previous global analyses 

In Ref. [25], the global analysis of J/ip suppression in photon-nucleus (E691) and 
hadron-nucleus (E537, NA3, E672, E772) reactions is carried out by Gerschel and 
Hiifner (GH). Excluding the NA3 measurements from the global fit, they found that 
^j/vN ^ i 0-3 mb. A somewhat larger but compatible result, (J = 7.3 ± 0.6 mb, 
was later determined from the NA38 and E772 measurements of J/ip production in 
proton-nucleus collisions by Kharzeev et al. (KLNS) in Ref. [26]. These results are 
significantly larger than our present estimate. Let us try to clarify the origin of this 
discrepancy. 

The E537, NA3 and E672 data were analysed in both the GH work and the present 
paper. However, although the values extracted by GH from the NA3 measurements 
are small (around 3 mb) and compatible with our result, these data were not included 
in the GH global fit. Moreover, the E691 photoproduction data, which gave (Tj^^^ — 
6.0 ± 1.6 mb in GH, were not selected in this work (sec Sect. 4). More importantly, 
the most accurate estimate of the J/ip-N cross section found in [25] is due to the 
E772 measurements {cr j^^^ = 6.25 ± 0.5 mb). However, as already discussed, these 
data suffer from a bias in the experimental acceptance which systematically increases 
the nuclear dependence of the J/t/j suppression. The corrected E772 results would 
then yield a somewhat smaller value for the J/i/j-'N cross section, which should be 
compatible with (Tj^^-^ = 2.7 ± 0.9 mb presently extracted from the more recent E866 
experiment in the same kinematic range. Last but not least, many new experimental 
data (chronologically NMC, NA50, E866, PHENIX, and HERA-B) came out since the 
analysis of Ref. [25]. Surprisingly, all these results lead to J/t/j-N cross sections smaller, 
or even much smaller (see Table 2), than the estimate quoted in GH. 

Comparing our analysis with that of KLNS, we can remark that the cTj^^^ currently 
extracted from NA38 data (cr^^^j^ = 5.5±0.7 mb) is compatible with the value estimated 
in [26], CTj^^-^ = 6.2 ± 0.3 mb. However, the strong disagreement between the data and 
the theoretical calculations (x^/ndf = 3.4) weakens the confidence in this estimate 
and marginally affects the result from the global fit because of the smaller weight; see 
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Eq. (5). 



7.2 Comparing with phenomenological estimates 

Let us now compare our result with phenomenological expectations. The J/ip-N cross 
section has been computed in a generalized vector dominance model (VDM) in [10]. 

It is found to be 2.8 < crj^^ (v*j/v,jv = ^0 GeV^ < 4.1 mb, with a smooth energy 
dependence 

where the exponent of the power-law behaviour A ~ 0.4 is determined from the energy 
dependence of the J/ip photoproduction cross section. Given that the bulk of the data 
points analysed here he in a range y/sJJ^ ~ 10-20 GeV (see Sect. 4.1), the VDM 
estimate in this energy window then reads 

2.8-3.7 < a™ [,/^ = 10-20 Gev) < 4.1-5.4 mb, (19) 

in perfect agreement with the above extracted value cr ^^^^^ — 3.4 ± 0.2 mb. Note 
that a similar energy dependence is expected in the pQCD approach by Bhanot and 
Peskin [7], where the exponent A is now directly related to the behaviour of the gluon 
distribution in the nucleon at small x. However, it appears that the magnitude of the 
J/ip-N cross section in this approach is roughly cTj^^^^ — 1.5-2 mb in the energy domain 
of interest [9], that is a factor of 2 smaller than what was found here. Finally, let us 
mention that the non-perturbative estimate cTj^^^ — 3.6 mb determined in Ref. [46] 
matches our result well. 

In comparing our result with theoretical expectations, however, we should bear in 
mind that J/ ip production may also come from the radiative decays of higher states such 
as Xc and ip' which, because of their larger size, are expected to interact more strongly 
with the nuclear medium than does the J/ip. In that respect, the measurement of the 
nuclear dependence of Xc production together with a more accurate determination of 
the Xc and ijj' contributions to J/ijj production, possibly performed by the NA60 and 
HERA-B experiments, is highly desirable (see [47] for a discussion on this issue). 



7.3 Energy dependence 

We have implicitly assumed so far that the J/ip^s are produced instantaneously in 
the nucleus. This picture is obviously too naive because of the typical proper time 
Tj^^ ~ (m^, — nij^^)^^ ^ 0.3 fm needed to form the J/ip state [48,49]. At high 
energy, the Lorentz-contracted nuclear radius, -R/7, becomes smaller than r^^^; what 
propagates through the nuclear medium is thus rather a compact cc state than a fully 
formed meson. Several models based on colour transparency have been proposed to 
properly describe the propagation of an expanding cc dipole in nuclei (see e.g. [49,50]). 
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In our approach - which we want to keep as model-independent as possible - as well as 
in Refs. [25,26], such formation time effects are not addressed. Hence, the constant cross 
section cTj^^^ that enters Eq. (5) has to be seen and understood rather as an effective 
absorption parameter (resulting from the average of the cc and charmonium interaction 
with nucleons) than the genuine J/ip-'H cross section [48], despite the above-mentioned 
good agreement between our extracted value and the theoretical estimates. 
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Figure 5: The J/ip-'H cross section extracted from each data set as a function of x^. 
The solid (dashed) vertical bars indicate the (rescaled) la error, 5a (^^j/^n)' 
the dashed band represents the value obtained from the global fit of all data. 



Since the J/ip-N energy, Eq. (17) (as well as the boost, 7 



m 



/2m 



j^Xj, in the 



we could - at least 
/(xj. In 



a 



nucleus rest frame), is directly related to the momentum fraction Xj, 
in principle - expect the cross section to be a scaling function of Xj , - ^^^i^ 
QCD, a compact cc dipole does not interact strongly with the nuclear matter because 
of colour transparency [48,51]. Consequently, /(xj should be a c?ecreasm^ function as 
X2 tends to (i.e. in the high energy limit); yet this trend could somehow be balanced 
by the energy dependence of Eq. (18), which predicts a^^^^ ~ ^377777^ ~ xl 



-A/2 



j/i>N 



In order to investigate the possible energy dependence of the cc-N interaction, 
we plot in Fig. 5 the values extracted from the various data samples as a function 
of X2. The discrepancy between the E537/NA38/NA50 and the NA3/E866/HERA-B 
experiments already mentioned in Section 5 is clearly visible around ^0.1. Besides 
this observation, no real Xj-dependence of aj^^^ can be inferred from Fig. 5, despite the 
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different dynamics when going from ~ 0.2 down to the ~ 10~^ fraction probed 
at RHIC^ (this is less true when the EKS nPDF is used, see Sect. 6.2). Needless 
to mention that having more precise d-Au data at RHIC would significantly help to 
clarify the picture. On a longer term, it would also be interesting to extract similarly 
J/V'N from the J/ip suppression in proton-nucleus collisions at even higher energy (i.e. 
LHC), so as to eventually signal the onset of colour transparency. Finally, the lack of 
x^ dependence in the data justifies a posteriori our choice to perform a global fit of the 
existing measurements with a constant cross section, irrespective of the energy of the 
reaction. 



8 Summary 

In summary, all available data on J/ ip hadroproduction and muoproduction on nuclear 
targets have been analysed within a Glauber model. These measurements allow for 
the extraction of the J/ip-'N inelastic cross section from each experiment. Fitting data 
altogether leads to cTj^^-^ = 3.4 ± 0.2 mb (x^/ndf = 1.5). A similar analysis, including 
nuclear corrections in the PDFs, yields comparable results (c°^^ = 3.5 ± 0.2 mb, 
X^j^g/ndf = 1.4). We stress however that a much stronger gluon shadowing (such as 
the one given in nDSg or EKS) leads to a somewhat larger value, even though such a 
significant antishadowing seems disfavoured by the E866 J/ t/j data. 

Remarkably, our result is significantly lower than previous similar estimates [25,26]. 
This is partly the consequence of the new data available (in particular the E866, HERA- 
B and, to a lesser extent, the PHENIX measurements) which have appeared since 
the analyses of Rcfs. [25,26]. The cross section found in this paper is rather close to 
phenomenological expectations for the interaction oia. J/ip state with nucleons. Finally, 
the present uncertainties of the PHENIX data do not allow a statistically significant 
energy dependence of (7j,^^ to be observed. 
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